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摘要 采用 ANSYS 有 限 元 法 建立 三 维 模型 ,研究 了 基体 热膨胀 系数 和 弹性 模 量 对 复合 材料 界面 径 向 、 轴 向 和 周 向 残余 应 力 
大 小 和 分 布 的 影响 。 结 果 表 明 , 界面 残余 应 力 的 改变 量 均 与 基体 热膨胀 系数 的 改变 量 成 正比 , 且 随 着 基体 热膨胀 系数 的 减 
小 界面 周 向 残余 应 力 沿 着 纤维 周 向 分 布 的 不 均匀 性 减弱 ; 界面 残余 应 力 变化 的 总 趋势 是 随 着 基体 弹性 模 量 的 增 大 而 增 大 ， 
但 是 增加 的 幅度 不 断 减 小 。 
关键 词 复合 材料 ， 热 残余 应 力 ， 有 限 元 法 ， 基 体 性 能 
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ABSTRACT The magnitude and distribution of the interfacial stresses along the radial-, axial- and cir- 
cumferential directions at fiber side of fiber reinforced composites SiC:/ Ti-6Al-4V with the change of ther- 
mal expansion coefficient and elasticity modulus of the matrix are analysed by three-dimensional model 
established by using finite element method. The results show that the change of interfacial residual stress 
is proportional to that of thermal expansion coefficient of matrix, and the no-uniformity of distribution of 
the residual stress along the circumferential direction of fiber decreases with the decrease of thermal ex- 
pansion coefficient. In addition, interfacial residual stress usually increases with the increase of elasticity 
modulus of matrix, but the increase of which reduces gradually. 

KEY WORDS composites, thermal residual stress, finite element method, matrix properties 


金属 基 复 合 材 料 (MMCs) 具 有 较 高 的 比 强 度 、 沿 轴 向 方向 的 残余 压 应 力 使 其 在 纵向 载荷 作用 下 


比 刚度 、 比 模 量 、 极 佳 的 蠕 变 和 疲劳 性 能 以 及 良好 ”的 承载 能 力 增强 环 向 的 残余 拉 应 力 使 得 复合 
的 高 温 性 能 , 在 航空 、 航 天 、 电 子 及 运输 等 高 技术 “材料 基体 中 产生 一 些 垂 直 于 界面 的 径 向 裂纹 , 特 


领域 有 1) 


一 二 


阔 的 应 用 前 景 吧 。 但 是 , 由 于 SiC 纤 维和  ” 别 是 当 纤维 间距 较 小 时 “。 因 此 , 全 面 了 解 复合 


Ti 合金 基体 热膨胀 系数 有 和 较 大 差异 , 当 SiCyTi 复 。 材料 热 残 余 应 力 的 各 种 影响 因素 、 热 残余 应 力 的 
合 材 料 从 制备 态 、 热 处 理 或 二 次 加 工 温度 冷却 到 ”状态 及 分 布 规律 , 对 复合 材料 强度 估算 和 寿命 预 


室温 时 在 材料 中 产生 热 残余 应 力 *% 影响 其 力学 。 测 等 有 重要 意义 。 
。 例 如 , 纤维 /基体 界面 的 径 向 残余 压 应 力 使 本 文 以 Ti-6Al-4V 合金 为 基础 , 以 最 常见 的 纤维 
其 在 横向 载荷 作用 下 界面 的 脱 粘 延迟 0 纤维 中 ， 四 方 和 六 方 排 布 的 代表 性 体 元 为 模型 (纤维 体积 分 


性 能 


数 35%), 研究 SiCyTi 基 复合 材料 从 高 温 制 备 态 冷却 


至 室温 时 , 基体 热膨胀 系数 和 弹性 模 量 的 变化 对 


* 国家 自然 科学 基金 51271147 资 助 项 目 。 
A pe ee 、 
2015 年 6 月 5 日 收 到 初稿: 2015 年 12 月 17 日 收 到 修改 稿 SiCyTi 复 合 材 料 界面 径 向 i 向 和 周 向 热 残 余 应 力 
本 文联 系 人 : 娄 菊 红 大 小 和 分 布 的 影响 , 并 分 析 这 些 热 残 余 应 力 与 复合 
DOI: 10.11901/1005.3093.2015.162 材料 宏观 响应 之 间 的 关系 。 
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504 材 料 而 


1 有 限 元 模拟 

所 用 的 材料 为 SiCyTi-6AL4V 复合 材料 , 设 基体 
为 理想 弹 一 塑性 应 变 硬 化 材料 , 服从 Von Mises 屈服 
准则 , 在 不 同 温 度 下 热膨胀 系数 和 弹性 模 量 的 原始 
值 与 变化 值 分 别 列 于 表 1 和 表 2"* 设 纤维 为 线 弹 
性 材料 , 直径 为 100 pm, 其 性 能 列 于 表 3""。 
复合 材料 中 常见 的 纤维 排 布 方式 有 四 方 排 布 
和 六 方 排 布 , 其 二 维 截面 示意 图 和 对 应 的 有 限 元 
模型 如 图 1 和 图 2 所 示 。 模 型 在 Z 方 向 的 高 度 为 
1000 um,“0 ”表示 XOY 坐标 平面 内 从 义 方向 到 Y 
方向 旋转 的 角度 ， 表示 基体 沿 着 界面 的 各 个 位 
置 ; “2Z” 表 示 沿 纤维 长 度 方向 的 尺寸 。 模 拟 过 程 中 
假定 纤维 和 基体 处 于 理想 粘 合 状态 , 其 边界 条 件 见 
文献 [15]。 


2 结果 和 分 析 
3 给 出 了 不 同 热膨胀 系数 下 纤维 一 侧 界面 残 
表 1 Ti-6A1-4V 基体 热 脱 


完 学 


报 30 卷 


余 应 力 沿 纤维 周 向 (纤维 四 方 排列 9=0°~45° 位 置 和 
纤维 六 方 排列 的 0=0°~30? 位 置 ) 在 Z=500 wm 平面 上 
的 分 布 。 可 以 看 出 , 无 论 是 四 方 还 是 六 方 排 布 , 界面 
径 向 、 轴 向 和 周 向 残余 应 力 的 改变 量 均 与 基体 热 脱 
胀 系数 的 改变 量 成 正比 。 在 0° 处 纤维 四 方 和 六 方 


1 三 种 纤维 排 布 方式 的 二 维 截面 示意 图 (单位 : mm) 


Fig.1 Two-dimensional cross-sections of different fiber ar- 


rays (unit: mm) (a) square array, (b) hexagonal array 


系数 的 原始 值 及 变化 值 (107C) 
Table 1 Original data and change value of thermal expansion coefficient for Ti-6Al-4V matrix (107C) 


TNC) A 0 A OQ 
20 8.8 7.04 3.28 3.352 
200 9.4 7.52 5.64 3.76 
400 10.3 8.24 6.18 4.12 
600 10.8 8.64 6.48 4.32 
800 J1S 9.2 0.9 4.6 


表 2 Ti-6AL4V 基体 弹性 模 量 的 原始 值 及 变化 值 (GPa) 
Table 2 Original data and change value of elastic modulus for Ti-6Al-4V matrix(GPa) 


ey a a a 
20 114.0 228.0 342 456 
200 103.8 207.6 311.4 415.2 
400 92.6 185.3 277.8 370.4 
600 76.4 152.8 229.2 305.6 
800 62.8 125.6 188.4 251.2 


表 3 国产 SiC 纤 维 的 性 能 
Table 3 Main properties of domestic SiC fiber 
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至 


图 2 不 同 纤维 排 布 方式 的 


了 限 元 模型 


Fig.2 Finite element model for different fiber arrays (a) 


square array, (b) hexagonal array 
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排 布 情况 下 基体 热膨胀 系数 每 减 小 0.2a, 界面 径 
向 、 残 余 压 应 力 分 别 减 小 85 MPa 和 74 MPa, 且 当 室 
温 下 基体 热膨胀 系数 小 于 纤维 热膨胀 系数 时 界面 径 
向 残余 应 力 变 为 拉 伸 状态 (图 3al 和 a2), 严重 破坏 纤 
维 /基体 界面 结合 , 甚至 产生 沿 界面 的 环 向 裂纹 ; 界 


于 


油 向 残余 压 应 力 分 别 减 小 201 MPa 和 200 MPa， 


而 且 与 径 向 残余 应 力 的 情况 相同 , 当 基体 热膨胀 系 
数 小 于 纤维 热膨胀 系数 时 界面 轴 向 残余 应 力也 变 为 
拉 伸 状态 (图 3bl 和 b2), 使 复合 材料 (纤维 ) 承 受 纵 问 


拉 伸 载荷 的 能 力 减 弱 ; 界面 周 向 残余 压 应 力 分 别 减 
小 51 MPa 和 66 MPa, 甚至 变 为 拉 伸 状态 (图 3cl 和 


c2), 这 


了 利于 防止 垂直 于 界面 的 径 向 裂纹 。 在 扩 


45°( 四 方 排 布 ) 或 9=30°( 六 方 排 布 ) 处 纤维 四 方 和 六 
方 排 布 情况 下 基体 热膨胀 系数 每 减 小 0.2a, 界面 径 
向 残余 压 应 力 分 别 减 小 61 MPa 和 66 MPa, 减 小 的 


0 [7 站 fas™ YY 一 一 一 一 v 一 一 v 一 v 一 | 
| ol 
. a Om A va A ee 生 让 
在 
全 -50 人 人 人 一 人 
三 -一 全 一 一 六 50 下 人 一 全 一 全 一 一 全 一 一 全 一 一 全 一 一 全 一 一 全 一 一 全 一 一 1 
9 0 
$B -100 一 一 一 -一 | 蕊 
2 Pd 蕊 -100 下 
时 .15 一。 一 。 BB ee 一 一 一 一 一 ?一 “| 
避 | & 
国 
可 -150 上 
200 一 
一 一 一 下 
250[” 200 | 一 = 一 一” 
到 | | | 
10 20 30 40 0 10 20 30 
0/(") 0/(") 
[oy™ -一 YY- 一 YY- 一 YY 一 YY 一 YY- 一 YY 一 | | 画 " 一 YY 一 YY 一 YY 一 YY 一 YY 一 YY 一 YY 一 | 
Ee Hao a 一 RE 全 一 民 一 及 CA 
全 -100 上 他 -100 
过 仁和 三 一 人 人 人 人 人 人 一- 全 一 一 4 
出 -200 上 | % -200 
$s vg 
加 咏 
总 -300F 再 -300 
< eo—e—e—e—e—e—e— ee—4 < | ee—e6—e0—e—e0—e—e—e—4 
-400 上 -400 
-500 上 -500 
一 人 
0 10 20 30 40 0 5 10 15 20 25 30 
0 
/0 9/() 
| | 
0 Hel) 0 上 
要 图 a, @ % 全 Ce A a 而 Dd % 全 一 人 更 化 
人 一 一 全 一 全 G 
a 全 一 一 全 [ey 
50 上 一 和 
三 人 三 -50 全 一 一 全 一 一 全 一 一 全 一 人 一 人 人 人 人 1 
只 ) 一 e@ 加 
0 -100 上 “ee—e 加 
加 ”人 号 -100 
3 | 一 ee By 
8 -150f Me sd 9 0—e—e——e ee ee | 
亚 a 工 
L! 上 
3 ~、 150 
| 
es nn 四 四 二 
-250 上 | ee _200 1 1 1 下 
0 10 20 30 40 0 5 10 15 20 25 30 
0/() 0/() 
图 3 不 同 基体 热膨胀 系数 下 纤维 一 侧 界面 残余 应 力 的 分 布 


Fig.3 Distributions of interfacial residual stresses in fiber with different coefficient of thermal expansion 
for matrix (a)radial residual stress (b) axial residual stress (c) hoop residual stress (“1”represents fi- 
ber square array and “2”represents fiber hexagonal array) 


一 个 一 


506 


系数 越 接 近 径 向 残余 应 力 沿 纤 维 周 向 的 不 均匀 程 
度 越 小 (图 3al 和 a2); 界面 轴 向 残余 压 应 力 分 别 减 
小 204 MPa 和 199 MPa 时 , 变化 幅度 与 0=0° 位 置 基 


aN 
~ 


材 料 研究 学 报 


幅度 比 9=0° 位 置 的 小 , 这 说 明基 体 与 纤维 的 热膨胀 


本 相等 ; 界面 周 向 残余 压 应 力 分 别 减 小 90 MPa 和 
71 MPa, 减 小 幅度 比 9=0° 位 置 的 大 , 即 随 着 基体 热 
膨胀 系数 的 减 小 周 向 残余 应 力 沿 纤维 周 向 的 不 均匀 


程度 减 小 


o 


4 给 出 了 不 同 弹性 模 量 下 纤维 一 侧 界 面 残余 
应 力 沿 纤维 周 向 (纤维 四 方 排 列 9=0°~45° 位 置 和 纤 


维 六 方 排列 的 0=0°~30° 位 置 ) 在 Z=500 wm 平面 上 的 


分 布 。 可 以 看 出 , 随 着 基体 弹性 模 量 的 增 大 (基体 变 
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30 卷 


趋势 是 径 向 、 轴 向 和 周 向 残余 压 应 力 均 增 大 , 而 且 
基体 弹性 模 量 从 开始 增加 第 一 个 1 倍 的 鼠 ( 即 


弹性 模 量 从 El 增 为 妃 ) 时 界面 


残余 应 力 变化 量 最 


大 , 增加 第 二 个 1 倍 的 豆 ( 即 红 


性 模 量 从 到 增 为 


也 ) 时 界面 残余 应 力 变化 量 次 之 , 增加 第 三 个 1 倍 


口 


的 吾 ( 即 弹性 模 量 从 已 增 为 劲 ) 时 界面 残余 应 力 
变化 量 最 小 。 对 于 纤维 四 方 排列 , 在 0=0° 和 6=45 


位 置 随 着 基体 弹性 模 量 的 不 断 增 加 , 界面 径 向 残余 
压 应 力 的 增 量 依次 为 168 MPa、6 MPa、-27 MPa 和 


149 MPa、141 MPa、85 MPa( 图 4al), 即 弹性 模 量 从 


且 增 为 及 时 径 向 残余 压 应 力 在 


0=0° 和 45° 位 置 都 大 


大 提高 , 但 是 随 着 弹性 模 量 的 进 


i te 本 a i, 
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应 力 仅 在 G45? 位置 有 所 提高 而 在 G0° 位 置 变 化 不 大 ， 
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图 4 不 同 基体 弹性 模 量 下 纤维 


侧 界面 残余 应 力 的 分 布 


Fig.4 Distributions of interfacial residual stresses in fiber with different elastic modulus for matrix (a) ra- 


dial residual stress (b) axial residual stress (c) hoop residual stress (“1” represents fiber square ar- 


ray and “2”represents fiber hexagonal array) 
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使 得 界面 结合 的 薄弱 环节 由 0=45° 位 置 逐 渐 转 向 
9=0° 位 置 ; 界面 轴 向 残余 压 应 力 的 增 量 依次 为 
282 MPa、 110 MPa、 67 MPa 和 307 MPa、 203 MPa、 
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越 接 近 热 残余 应 力 沿 纤 维 周 向 分 布 的 均匀 性 越 好 。 
2. SiCyT 复合 材料 纤维 一 侧 界 面 径 向 、 轴 向 和 
周 向 热 残 余 压 应 力 变 化 的 总 趋势 是 随 着 基体 弹性 模 


88 MPa( 图 4b1), 即 界面 轴 向 残余 压 应 力 始终 随 着 基 
体 弹 性 模 量 的 增加 而 增加 , 但 是 当 弹 性 模 量 增 量 恒 
定 (1 倍 的 囊 ) 时 界面 轴 向 残余 压 应 力 的 增幅 逐渐 减 
小 , 且 沿 纤维 周 向 的 不 均匀 程度 加 大 ; 界面 周 向 残余 
压 应 力 的 增 量 依次 为 92 MPa、-4 MPa、40 MPa 和 
211 MPa、230 MPa、13 MPa( 图 4c1), 可 见 基体 弹性 
模 量 的 增加 使 界面 周 向 残余 压 应 力 沿 纤维 周 向 的 不 
均匀 程度 迅速 提高 ( 周 向 残余 应 力 在 0=0° 位 置 远 远 
小 于 0=45° 位 置 ), 即 高 的 基体 弹性 模 量 很 容易 引起 
945° 位 置 的 径 向 裂纹 。 对 于 纤维 六 方 排列 , 在 0° 
和 0=30° 位 置 随 着 基体 弹性 模 量 的 不 断 增 加 , 界面 径 
向 残余 压 应 力 的 增 量 依次 为 173 MPa、72 MPa、6 MPa 
和 164 MPa、113 MPa、49 MPa( 图 4a2), 与 纤维 四 方 
排 布 情况 下 界面 径 向 残余 压 应 力 的 变化 类 似 , 当 基 
体 弹 性 模 量 超过 甩 后 界面 结合 的 薄弱 环节 由 人 30? 
位 置 逐 渐 转 向 0° 位 置 , 但 是 径 向 残余 压 应 力 沿 纤 
维 周 向 的 分 布 相对 均匀 ; 界面 轴 疝 残余 压 应 力 的 增 
量 依次 为 300 MPa、144 MPa、77 MPa 和 301 MPa、 
165 MPa、73 MPa( 图 4b2), 即 随 着 基体 弹性 模 量 的 增 
加 G0=0°* 和 30° 位 置 界面 轴 向 热 残 余 压 应 力 的 增幅 基 
本 一 致 , 说 明 轴 向 热 残余 压 应 力 沿 纤 维 周 向 的 均匀 
分 布 不 因 基 体 弹 性 模 量 的 变化 而 变化 , 且 弹 性 模 量 
每 次 的 增 量 恒定 (1 倍 的 Bl) 时 界面 轴 向 残余 压 应 力 增 
幅 为 上 一 次 增幅 的 1/2 左右; 界面 周 向 残余 压 应 力 的 
增 量 依次 为 159 MPa、53 MPa、69 MPa 和 172 MPa、 
95 MPa、9 MPa( 图 4c2), 可 见 基 体 弹 性 模 量 的 增加 使 
界面 周 向 热 残 余 压 应 力 沿 纤维 周 向 的 不 均匀 程度 先 
增加 后 减 小 , 这 意味 着 弹性 模 量 的 进一步 增加 (高 于 
BE4) 可 能 导致 界 面 周 癌 热 残余 压 应 力 的 最 大 值 由 全 
30° 位 置 变 为 0=0° 人 位置。 可见, 随 着 基体 弹性 模 量 的 
增加 纤维 六 方 排 布 时 周 向 残余 应 力 的 变化 比 四 方 排 
布 时 平缓 得 多 , 且 六 方 排 布 时 周 向 热 残余 压 应 力 峰 
值 越 来 越 小 于 四 方 排 布 时 的 应 力 峰 值 。 总 之 , 基体 
弹性 模 量 的 提高 , 一 方面 使 得 径 疝 与 轴 疝 热 残 余 压 
应 力 增 大 有 利于 纤维 /基体 界面 结合 和 纵向 拉 伸 性 
能 , 另 一 方面 使 得 周 向 热 残余 压 应 力 增 大 可 能 造成 
基体 更 大 的 塑性 流 变 甚至 径 向 裂纹 "5。 


3 结 论 


1.SiCyT 复合 材料 纤维 一 侧 的 界面 径 向 、 轴 向 
和 周 向 残余 压 应 力 都 随 着 基体 热膨胀 系数 的 减 小 而 
减 小 , 甚至 变 为 拉 伸 状态 , 且 基体 与 纤维 热膨胀 系数 


< 


量 的 增 大 而 增 大 , 但 是 随 着 基体 弹性 模 量 的 增加 热 

残余 压 应 力 增加 的 幅度 在 减 小 。 
3. 随 着 SiCyT 复合 材料 基体 弹性 模 量 的 增加 纤 

维 六 方 排 布 时 热 残 余 应 力 的 变化 比 四 方 排 布 时 更 为 

平缓 , 且 六 方 排 布 时 界面 周 向 热 残 余 压 应 力 峰 值 越 

来 越 小 于 四 方 排 布 时 的 应 力 峰 值 。 

参考 文献 


1 YUAN Meini, YANG Yanqing, LUO Xian, ZHANG Rongjun, Mi- 


crostress distribution of titanium matrix composites, Chinese Jour- 
nal of material research, 22(4), 389(2008) 

( 原 梅 妮 , 杨 延 清 , 罗 贤 , 张 荣军 , 钛 基 复 合 材料 中 的 微 区 应 力 分 
布 , 材料 研究 学 报 , 22(4), 389(2008)) 


2 D.B. Miracle, Metal matrix composites - From science to techno- 


logical significance, Composites Science and Technology, 65, 2526 
(2005) 

3 Y.Q. Yang,Y. Zhu,Y. Chen, Q. G. Zhang, J.M. Zhang, Processing 
and property of SiC fiber reinforced Ti-matrix composite, Rare Met- 
al Materials and Engineering, 31, 201(2002) 

4 HU Hengzhi, Challenge and development trends to future aero-en- 
gine materials, Journal of aeronautical materials, 18(4), 52(1998) 
( 传 恒 志 , 未 来 航空 发 动机 材料 面临 的 挑战 与 发 展 趋向 , 航空 材料 
学 报 , 18(4), 52(1998) ) 

5 M.Kuntz,B. Meier, G. Grathwohl, Residual Stresses in Fiber-Rein- 


forced Ceramics due to Thermal Expansion Mismatch, Journal of 
the American Ceramic Society, 76(10), 2607(1993) 

6 MaZhi-jun, Yang Yan-qing, Zhu Yan, Chen Yan, Effect of fiber dis- 
tribution on residual thermal stress in titanium matrix composite, 
Trans. Nonferrous Met. Soc. China, 14(2), 330(2004) 

7 S. W. Warrier, P Rangaswamy, M. A. M. Bourke, S. Krishnamur- 
thy, Assessment of the fiber/matrix interface bond strength in SiC/ 
Ti-6Al-4V composites, Materials Science and Engineering A, 259, 
220(1999) 

8 J.A.Sherwood, H.M. Quimby, Micromechanical modeling of dam- 
age growth in titanium based metal-matrix composites, Computers 
& Structures, 56(2-3), 505(1995) 

9 M.P. Thomas, M. R. Winstone, Longitudinal yielding behavior of 
SiC-fibre-reinforced titanium-matrix composites , Composites Sci- 
ence and Technology, 59, 297(1999) 

10 A. Brunet, R. Valle, A. Vassel, Intermetallic TiAl-based matrix com- 

posites: investigation of the chemical and mechanical compatibility 

of a protective coating adapted to an alumina fibre, Acta Materialia, 

48, 4763(2000) 

HUANG Bin, Study on the interface, microstructure and property 

of SICYTi composites, Dissertation for the degree of PhD , North- 

western Polytechnical University(2010) 

( 黄 斌 , SiCVTi 基 复合 材料 的 界面 及 组 织 性 能 研究 , 博士 学 位 论 

文 , 西北 工业 大 学 (2010)) 

12 S. Schuler, B. Derby, M. Wood, Matrix flow and densification dur- 


ing the consolidation of matrix coated fibers, Acta Materialia, 48, 


一 
fk 


CN 


一 个 一 


0 


202303.10742v1 


chinaXiv 


~- 


508 材 料 研 
1247(2000) 
13 Z.J. Ma Y Q. Yang, X. H. Lu, X. Luo, Y. Chen, The effect of ma- 


14 


1s 


trix creep property on the consolidation process of SiC/Ti-6Al-4V 
composite, Materials Science and Engineering A, 433, 343(2006) 

R. J. Zhang, Y. Q. Yang, W. T. Shen, Preparation and tensile test of 
SiC fiber fabricated by three-stage chemical vapour deposition, J In- 
org Mater, 25(8), 840(2010) 

LOU Juhong, YANG Yanqing, LUO Xian, CHEN Yan, Finite ele- 


ment analysis for effects of fiber arrays on the residual stress of 


ChinaXiv 合 作 期 刊 


学 报 30 卷 


SiC fiber reinforced Ti-matrix composite, Rare metal materials and 
engineering, 40(2), 243-246(2011) 

( 娄 菊 红 , 杨 延 清 , 罗 贤 , 陈 彦 ,纤维 排 布 方式 对 SiC/Ti 材料 残余 
应 力 影响 的 有 限 元 分 析 , 稀有 金属 材料 与 工程 , 40(2), 243(2011)) 


16 LUO Xian, Study on the fabrication and properties of SiC fiber re- 


inforced copper-matrix composites, Dissertation for the degree of 
PhD, Northwestern Polytechnical University(2008) 

( 罗 贤 ,SiC 纤 维 增强 Cu 基 复 合 材料 的 制备 及 性 能 研究 , 博 了 
位 论文 , 西北 工业 大 学 (2008)) 


由 


